Zero modes and charged Skyrmions in graphene bilayer 
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We show that the electric charge of the Skyrmion in the vector order parameters that characterize 
the quantum anomalous spin Hall state and the layer-antiferromagnet in a graphene bilayer is four 
and zero, respectively. The result is based on the demonstration that a vortex configuration in two 
broken symmetry states in bilayer graphene with the quadratic band crossing has the number of zero 
modes doubled relative to the single layer. The doubling can be understood as a result of Kramers' 
theorem implied by the "pseudo time reversal" symmetry of the vortex Hamiltonian. Disordering 
the quantum anomalous spin Hall state by Skyrmion condensation should produce a superconductor 
of an elementary charge 4e. 



Bernal-stacked graphene bilayer is an interesting ex- 
ample of a two-dimensional system of electrons with the 
quadratic band crossing. When the Fermi energy is tuned 
to the band crossing point, a simple argument shows 
that the non-interacting ground state is unstable, and 
the gap in the spectrum opens at an infinitesimal repul- 
sive interaction. 1] Indeed, there have been several re- 
cent experiments that find a broken-symmetry insulat- 
ing ground state and the concomitant spectral gap in 
graphene bilayers.jH-Q The exact nature of the insulat- 
ing ground state is presently a point of contention, since, 
at least for the Coulomb repulsion between electrons, the 
ground states with different broken symmetries may lay 
rather close in energy. 

[Ml 

Since the different broken symmetry states in bilayer 
graphene appear to be nearly degenerate, their topologi- 
cally nontrivial spatial configurations, such as vortices, 
merons, and skyrmions, which would represent stable 
thermal or quantum fluctuations, become highly relevant. 
These topological defects in single layer of graphene are 
known to carry non-trivial quantum numbers. [Tsj since 
the corresponding Dirac Hamiltonian is linear in momen- 
tum and closely related to the example introduced by 
Jackiw and Rossi in their classic study of charge frac- 
tionalization in two dimensions, [l^ In the single layer 
graphene, however, at least without a magnetic field, any 
breaking of symmetry requires an interaction strength 
which appears to be too high.[l7| In the bilayer, on the 
other hand, the Hamiltonian is quadratic in momentum 
and without the Lorentz symmetry. As a result, its spec- 
trum in the presence of topologically nontrivial configu- 
rations of the order parameters (masses) is unknown. In 
this note we discuss some global features of this spectrum 
and point out some of its physical consequences. 

We first consider the Hamiltonian for the spinless elec- 
trons in the bilayer graphene in the presence of a unit vor- 
tex in two, out of three, mutually anticommuting mass 
terms. These masses represent: 1) the state with a bro- 
ken layer- inversion symmetry (BLIS), and 2) the two val- 
ley ferromagnets (VFMi, VFM2) which break both the 
layer inversion and the time reversal symmetries. The 
forth possible mass, which commutes with the above 
three, is the quantum anomalous Hall (QAH) state. This 



is the basic Hamiltonian onto which the other problems 
can, as we show, be mapped. This Hamiltonian, while 
violating the physical time reversal symmetry, exhibits 
another hidden antiunitary symmetry, which we name 
"pseudo time reversal" . Unlike the true time reversal 
for spinless particles, the pseudo time reversal symme- 
try squares into — 1 and implies double degeneracy of the 
entire spectrum. When combined with the chiral symme- 
try, this further implies that the zero modes, if present in 
the spectrum, must come in Kramers pairs of equal chi- 
rality. We obtain the two zero modes analytically for a 
representative Hamiltonian with the amplitude of the or- 
der parameter increasing linearly with the radius. Once 
found this way, the zero modes cannot shift their en- 
ergy upon a deformation of the amplitude, as long as the 
Hamiltonian respects the chiral symmetry. Filling the 
doubly degenerate zero level fully thus leads to a unit 
electric charge e being bound to the vortex. 

The above result is then used to show that the 
Skyrmion texture in the vector order parameter that 
characterizes the quantum anomalous spin Hall (QSH) 
state in the bilayer graphene carries the electric charge of 
4e. In contrast, the charge of the Skyrmion in the vector 
order parameter that describes the layered antiferromag- 
net (LAF) is zero. Quantum disordering the QSH state 
by the condensation of Skyrmions could thus lead to an 
exotic superconducting state of four-electron composites, 
which seems to lie outside the usual BCS paradigm. In- 
stead, one may think of it as a higher dimensional analog 
of Frohlich mechanism for ideal conductivity in one di- 
mension. [18] 

Pseudo time reversal symmetry. - Assuming that it is 
the A-sublattices of the two graphene layers that arc cou- 
pled by the interlayer hopping t± , the low-energy Hamil- 
tonian for the bilayer is 



-ffo = Vi7i + V272, 



(1) 



with Vi — piP2/m, V2 = ipi — p'^)/2m, with pi as 
the momentum operators measured from the band cross- 
ing point, 2m = t^/t^, with t as the nearest neigh- 
bor intralayer hopping integral. [19j We chose the rep- 
resentation in which 71 = tJa (g) 0-2, 72 = ctq (8> ci .|20l| 
Ho acts on the four-component Dirac spinor = 



2 



{bi+,b2+,bi-,b2-)'^ , where the first index labels the 
layer, the second labels the band crossing points ±7?, 
and all the components are on the B-sublattices of the 
two layers. [iJl We have assumed spinless fermions, and 
will include the spin shortly. The above Hamiltonian 
can also be understood in the context of the single layer 
graphene, as the second-order term in the expansion at 
the Dirac point in powers of the momenta. Its form is 
dictated by the Cs^ symmetry of the honeycomb lattice. 

Let us first review some of the familiar symmetries 
of iifn.|20i] In the continuum limit, the translational in- 
variance is an exact (continuous) symmetry. It is gener- 
ated by the momentum-like operator 735 = (T3 (8)cro- The 
time reversal, which consists of the exchange of the two 
Dirac points and the complex conjugation K , is also an 
exact (discrete) symmetry. If we define 70 = ao ® (T3, 
73 = (Ti cr2, and 75 = cr2 (E) (72, then the time reversal 
operator may be written explicitly as It — ijij^K , and 
735 = *7375- Note that if = +1, since the real spin has 
not been included yet. Finally, 72 represents the opera- 
tion of exchange of the layers, and when accompanied by 
the inversion of one of the two axis, it is also an exact 
(discrete) symmetry of Hq. 

There are only four mass matrices that anticommute 
with Hq and would gap out the spectrum: 70, 73, 75, 
and i7i72. They represent, BLIS state, VFMi, VFM2, 
and the QAH state, respectively, quite similarly to the 
nomenclature in the single layer graphene. [l^ Only the 
first and the last state respect translational invariance, 
and all but the first state violate the time reversal invari- 
ance. (See Table 1.) This means that, in contrast to the 
single layer graphene, |2l| it is impossible to form a four- 
component bilayer Hamiltonian which would contain a 
vortex configuration in two of the masses and would si- 
multaneously obey the time reversal symmetry. 

We now introduce the notion of the pseudo time rever- 
sal invariance. Let us define the antilinear operator 



A = 735^4 



(2) 



The only mass term that fails to commute with A is the 
one for QAH state, whereas the other three, together with 
the Hamiltonian Hq, commute with it. (See Table 1.) We 
therefore choose an arbitrary pair of the first three masses 
to form a vortex, and study the electronic spectrum in 
its presence. The vortex Hamiltonian is defined as 



Hy = II0 + A(r)(73 cos -I- 75 sin < 



(3) 



where (r, (p) are the polar coordinates. While Hy now 
violates both the translational invariance (735) and the 
time reversal symmetry (It) separately, it respects the 
pseudo time reversal symmetry: A] = 0. Moreover, 
since the matrix 735 can be thought of as the generator 
of translations, it itself is odd under time reversal. This 
implies, however, that 



A' 



-1. 



(4) 



The existence of such an antiunitary symmetry is, by 
Kramers' theorem. [2^] sufficient for the spectrum of Hy 



State 


mass 


735 


It 


A 


BLIS 


70 


-1- 
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VFMi 


73 
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VFM2 


75 






+ 


QAH 


47172 


+ 







TABLE I: The properties of the broken layer-inversion sym- 
metry (BLIS), two valley ferromagnets (VFMi, VFM2), and 
the quantum anomalous Hall state (QAH) of the spinless 
fermions in graphene bilayer under translations (735), time 
reversal (It), and the pseudo time reversal symmetry (A). 



to be at least doubly degenerate. We emphasize that 
this degeneracy has nothing to do with the electron spin, 
which is not included into Hy, nor with the true time 
reversal symmetry, which is a) broken by the mass terms 
in Hy, and b) even if it were not, it squares to -1-1, and 
by itself therefore does not imply any degeneracy. 

The uncovered Kramers' degeneracy of the spectrum 
of Hy is independent of the particular representation. All 
four-dimensional representations of the four matrices in 
Hy are equivalent. In particular, one can always choose 



7i to be real for i = 1, 2, and imaginary for i = 3,5, |15 | 
so that the pseudo time-reversal operator becomes easy 
to discern: A = 7375/^. 

Zero modes. - The spectrum of Hy contains two states 
with exactly zero energy. Assume that near the origin 
the amplitude vanishes linearly, A(r) = cr + 0{r^), so 
that the vortex configuration has a finite core energy. 
Discarding the higher order terms in this expansion, we 
can define the linearized vortex Hamiltonian 



Hyjin ^ Ho + c(a;i73 + 2:275), 



(5) 



so that Hy = Hy^iin + 0{r^). Both Hy and its linearized 
version Hy^un obey the pseudo time reversal symmetry, 
as well as the chiral symmetry: 



{iJt,,7o} = {Hy^iin,-fa} = 0. 



(6) 



Since 70 is even under pseudo time reversal, if the spec- 
trum of Hy^iin contains a Kramers' doublet at zero, both 
states will be the eigenstates of 70 with the same eigen- 
value. One may then think of the original Hamilto- 
nian Hy as a continuous amplitude deformation of Hy^un 
which preserves the chiral and the pseudo time reversal 
symmetries, so that the zero modes cannot neither move 
nor split in the process. 

It therefore suffices to show that Hy^un has two zero 
modes. This is easiest to demonstrate in the momentum 
representation, in which Xi = id/dpi, and 



Hy 



d d 

«^73+«T^75+P^(72 cos29p+-fi sin26'p) (7) 
dpi dp2 



where {p. Op) are the polar coordinates in the momen- 
tum space. In writing the last expression we have also 
introduced the dimensionless variables, so that energy 
is measured in units of {c^ /2mY^^ , and length in units 
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of (2toc)~^/^. The linearized Hamiltonian can now be 
viewed as describing a double vortex in momentum space, 
and on rather general grounds one expects it to have two 
zero modes. [16| Indeed, they can be analytically com- 
puted: 

^i = {0,e^'''f{p),0,9{p)f, (8) 
^2^{0,g{p),0,e-''''f{p)f, (9) 

with 

gip)=^fp^Ky3ipy3), (10) 

with K2/3{z) as the modified Bessel function, A/" as a 
normalization factor, and f{p) = —g'{p)/p^- Both g{p) 
and f{p) decrease exponentially at large momenta, and 
both are regular at small momenta. We also confirm that 

*2 = A*i, (11) 

by recalling that under time-reversal. Op O^p = 6p + t: 
as well. Both zero-modes have negative chirality: 7o^i = 
—"^i- For an antivortex their chirality would be positive. 

Irreducible Skyrmion Hamiltonian. - The existence of 
two zero modes implies that when both states are oc- 
cupied the electrical charge of the vortex is unity. This 
has interesting consequences for the charge of Skyrmion 
excitations. Before finally restoring spin and considering 
a spin-Skyrmion, let us define an auxiliary irreducible 
Hamiltonian 

H{rn) = Ho + mi(r)73 + "12(^75 + ™3(^7o, (12) 

with the three-component vector rn{r) describing a 
Skyrmion configuration: mapping of the plane onto the 
two dimensional sphere S2, at which the vector rn{f) with 
a fixed unit magnitude lives, with the Pontryagin index 

P = / cPfrn{'r) ■ {dxTn{f} x dyrn{r)), (13) 
47r J 

equal to unity. It is straightforward to see, on the other 
hand, that the Pontryagin index can be rewritten as 

P = ^iJ2nirv)coseir,) - cose{R)Y,nir,)), (14) 
where, 

n{r^) = jdl- \7(j){f) (15) 

is the vorticity at the point f\,, i? is the boundary of the 
system, and to = (sin cos sin sin 0, cos 0) . Skyrmion 
is usually depicted as a vortex with n = 1 at the ori- 
gin, and with the angle 6 interpolating from 6(0) = to 
0{R) — TT. But, Eq. (14) implies that Skyrmion is topo- 
logically equivalent to a pair consisting of a meron with 
n = 1 and = at the origin, and of an antimeron with 
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FIG. 1: The elementary Skyrmion in the mass m, as a meron- 
antimeron pair. Both the meron and the antimeron carry a 
unit electric charge of the same sign. 

n = — 1 and = tt at some other point in space (Figure 
1). If we assume the meron and the antimeron to be spa- 
tially well separated, simple counting leads to the elemen- 
tary skyrmion's total charge to be 2e. Explicitly, a meron 
at the origin is topologically equivalent to the configura- 
tion TO = (xi, X2, —to), and an antimeron at some distant 
position L is equivalent to to = (xi— Li, — (a;2— .^^2), -Vrri). 
Since the zero modes of a vortex and an antivortex have 
opposite chiralities, the charges of such a meron and an 
antimeron are the same and equal to e. The charge of 
the elementary Skyrmion configuration is therefore 2e. 

Charge of spin- Sky rmions. - All is ready now to restore 
the electron spin. Let us consider two prominent candi- 
dates for the ground state of the bilayer graphene: QSH 
state, and LAF. Their representative masses are given by 
N ■ a Z7i72 and ./V ■ it €5 70, with N as the vector order 
parameter, [23 respectively. Both states spontaneously 
break the rotational symmetry, and the latter also breaks 
the physical time reversal symmetry. Consider then a 
Skyrmion configuration for the order parameter N{r): 

H{N) ^cTo(g)Ho + Nir) ■a®M, (16) 

with M — 47172, or M = 7o. The five eight-dimensional 
matrices ctq 71, <Ja ® I2 and a ® M featured in H{N) 
all mutually anticommute and square to unity, and thus 
form a representation of the Clifford algebra C(5, 0). [lH| 
C(5, 0), on the other hand, has only two inequivalent irre- 
ducible representations, which are both four dimensional, 
and differ in the sign of an odd number of matrices. This 
guarantees that the spin-Skyrmion Hamiltonian can be 
decomposed by a unitary transformation as follows: 

U^H{N)U ^ H{N)® H{N') (17) 

where N' = N or N' = —N. The plus (minus) sign 
corresponds to two equivalent (inequivalent) representa- 
tions of C(5, 0) appearing in the two blocks. The electri- 
cal charge of the spin-Skyrmion depends on the result of 
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this decomposition: the plus sign imphes that the spin- 
Skyrmion is equivalent to the sum of the two Skyrmions 
in each irreducible block, so that the electrical charge of 
the two adds up to 4e, whereas the negative sign implies 
the sum of a Skyrmion and an anti-Skyrmion, which leads 
to the charge zero. 

To determine the charge of the spin-Skyrmion we note 
that the trace of the product of the five matrices in H{N) 

Q = ^Tr Y[ ao (E> J^ Y[ ak®M (18) 

1=2,1 fe=l,2,3 

is invariant under a unitary transformation, and (5 = 4 
for M = 17172 and Q = for M = 70. Since the same 
trace equals four for N = N' and zero for = — A^', we 
find that the electrical charge is finite only in the QSH 
state, when it is 4e. 

Topological superconductivity. - Since the charged 
spin- Skyrmions should be bosons, inducing their conden- 
sation by tuning an appropriate interaction term in the 
Hamiltonian which reduces their cost in energy would 
produce a superconductor with the elementary charge of 
4e. Another possibility is doping the system with a suf- 
ficiently large density of charge carriers. The ensuing 
superconducting state cannot naturally be thought of in 
terms of Cooper pairing mediated by some attractive in- 
teraction. Rather, it is the result of topologically non- 
trivial configuration in the QSH order parameter being 
created in order to fill their zero modes with extra elec- 
trons and then slide through the system, as discussed by 
Frohlich in the pre-BCS era. Frohlich's mechanism for 
ideal conductivity in one dimension has been extended 
to and greatly elaborated in higher dimensions by Wieg- 
mann [1^, and discussed in the context of cuprates [23 
and, recently, chalogenidesjl^. The novel point in our 



example is that due to the quadratic band crossing the 
charge of the Skyrmion becomes doubled, which makes 
the topological mechanism for superconductivity more 
efficient, and the contrast with the BCS paradigm even 
more striking. 

Conclusion and summary. - We demonstrated that the 
charge of Skyrmion textures in a graphene bilayer with 
quadratic band crossing is doubled relative to its single 
layer equivalent. The basis for our result is the explicit 
derivation of two zero modes for the basic vortex Hamil- 
tonian with quadratic dispersion, and the identification 
of the hidden pseudo time reversal symmetry which pro- 
tects their degeneracy. It was argued that the condensa- 
tion of the Skyrmion textures of the QSH state in bilayer 
graphene upon doping leads to a non-BCS superconduc- 
tor with the flux quantized in units hc/Ae. 

The present derivation of the charge of the Skyrmion 
in the QSH state can be applied equally well to the sin- 
gle layer graphene Hamiltonian, in which case it yields 
the charge of 2e, in agreement with Our approach 
circumvents the use of the gradient expansion and of the 
concomitant Wess-Zumino-Witten term [g^], the applica- 
bility of which to systems with the quadratic band cross- 
ing does not seem obvious '285, but may be possible. [2^ 

The nature of the doubling discussed here is different 
from the one for the bilayer graphene in the magnetic 
field, [1^ or in one-dimensional systems with quadratic 
band crossing p8l |. when only the number of zero modes, 
and not the whole spectrum, is doubled, and the pseudo 
time reversal symmetry introduced here is absent. 
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